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The structural phase transitions and the fundamental band gaps of Mga;Zni_a;0 alloys are inves- 
tigated by detailed first-principles calculations in the entire range of Mg concentrations x, applying 
a multiple-scattering theoretical approach (Korringa-Kohn-Rostoker method). Disordered alloys 
are treated within the coherent potential approximation (CPA). The calculations for various crystal 
phases have given rise to a phase diagram in good agreement with experiments and other theoretical 
approaches. The phase transition from the wurtzite to the rock-salt structure is predicted at the Mg 
concentration of a; = 0.33, which is close to the experimental value of 0.33 — 0.40. The size of the 
fundamental band gap, typically underestimated by the local density approximation, is considerably 
improved by the self-interaction correction. The increase of the gap upon alloying ZnO with Mg 
corroborates experimental trends. Our findings are relevant for applications in optical, electrical, 
and in particular in magnetoelectric devices. 

PACS numbers: 61.50.Ks,81.30.Hd 
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I. INTRODUCTION 

In recent years, much effort has been devoted to re- 
search on ZnO; inspired mostly by its attractive proper- 
ties for optoelectronic applications This interest arises 
from specific properties, e. g. a large piezoelectric coeffi- 
cient, photoconductivity, and transparency in the visible 
and infrared wavelength regimes. The range of applica- 
tions of this semiconductor can be considerably extended 
by alloying. Prominent dopants are Co and in particu- 
lar Mg, on which we focus in this work. An increase 
of the Mg concentration can transform the crystal lattice 
from the wurtzite structure (WZ) of ZnO to the rock-salt 
structure (RS) of MgO. Accompanied by this structural 
phase transition is a substantial increase of the funda- 
mental band gap. The latter can be tuned from 3.35 eV 
to 7.7eVjii^ In view of a possible band-gap engineering, 
Mga;Zni_2,0 alloys may also be considered as suitable 
insulating spacers in magnetoelectronic devices, in par- 
ticular in magnetic tunnel junctions. 

According to the equilibrium phase diagramj^i^ the 
solid solution of Mg2,Zni_a;0 is of eutectic type at nor- 
mal conditions. It is characterized by an extensive sol- 
ubility of zincite in MgO (up to 33 mol%) and by a re- 
stricted solubility of MgO in ZnO (4mol%). The solu- 
bility depends strongly on experimental conditions and 
can be considerably increased at high temperatures and 
high pressuresi^i^ Non-equilibrium growth processes, like 
pulsed laser deposition (PLD)^i^'^ i^°i^^ and molecular 
beam epitaxial method o^^i^^i^^'^^ , allow to grow high- 
quality Mga;Zni_.^0 thin films for a large range of con- 
centrations X. For RS-Mg2;Zni_3;0, produced by PLD, 



a maximum solubility has been reported for x = 0.5^^ 
When increasing the Zn concentration, the RS and the 
WZ phases separate, and Mga;Zni_3;0 exhibits a WZ 
structure for x < OAJ^ The solubility limit and the 
phase formation in Mg^Zni^xO are strongly influenced 
by experimental conditions and by the substrate on which 
the Mga;Zni_2;0 film is grown. In general, alloying of 
ZnO and MgO proceeds by substituting Mg atoms by Zn 
atoms in the cubic RS structure and vice versa in the 
hexagonal WZ structure. 

The composition and the crystalline structure affect di- 
rectly the electronic properties of Mg^^Zni-^jO. Numer- 
ous absorption and photoluminescence spectroscopy ex- 
periments show that the fundamental band gap depends 
differently on the Zn concentration in the RS and WZ 
phaseSi ^^i^^ The width of the band gap increases most 
linearly with Mg concentration for both the RS and the 
WZ phase, but the slope in cubic Mg3;Zni_2,0 is about 
twice as large as in the hexagonal structure This ex- 
perimental finding clearly indicates that Mga;Zni_3;0 is 
a promising candidate for band-gap engineering. For 
instance in a magnetic tunnel junction, electrons that 
arc transmitted from one electrode to the other have to 
pass the nonconducting spacer The transmission prob- 
ability decays with spacer thickness and with the width 
of the spacer's fundamental band gapj ^^'^"^ As a conse- 
quence, the spin-dependent conductance, i. e. the tunnel 
magnetoresistance, could be tuned by varying the fun- 
damental band gap. Further, the magnetoresistance de- 
pends essentially on the properties of the ferromagnet- 
insulator interface, as was shown for Fe/MgO/Fe tunnel 
junctionSi ^^i^^ Hence, detailed knowledge of its geometric 
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structure is necessary, and our theoretical investigation 
of the bulk structural phases can be regarded as one step 
towards that goal. 

In addition to the extensive experimental work, there 
exist many detailed theoretical studies of both ZnO and 
MgO. Among them are Refs. [11123128,29,30,31,32,33.33, 
l35ll36l . using a variety of ab initio computational methods. 
There are, however, not many studies of Mgx'Zini-xO 
alloys, possibly due to the complex interplay of their 
electronic and geometric structures. Recently, thermo- 
dynamical stability and ordering tendencies of the al- 
loys have been carefully investigated using the cluster- 
expansion methodi^ Based on the parameterization of 
total energies for various alloy configurations, the lat- 
ter allows to study accurately structural properties, in- 
cluding short-range order (SRO) effectsi^i^ It is found 
that the RS ^ WZ transition occurs at x w 0.33, which 
is consistent with experiment (x = 0.33, Ref. |7|; simi- 
lar results were obtained by a pseudopotential method 
using supercells^) . The cluster expansion method was 
also used by Seko et al. for investigating phase tran- 
sitions, including vibrational effects through lattice dy- 
namics calculations 4i The authors demonstrate that the 
transition pressure decreases with increasing Mg content, 
which is explained as follows. Below the solubility limit 
of MgO in ZnO, the RS phase is energetically preferred 
to the WZ phase in MgO. Above the solubility limit, the 
configurational entropy increases by the transition from 
a mixed WZ-ZnO/RS-MgO to a RS single phase. 

A brief critical review of previous theoretical ab ini- 
tio investigations has shown that: (i) In ZnO, the en- 
ergy levels of the localized Zn-3d electrons arc found rel- 
atively high and thus close to the valence bands, result- 
ing in strong hybridization with 0-2p states. Since these 
hybridization effects have to be correctly taken into ac- 
count, the all-electron methods are preferred to the pseu- 
dopotential methods in which localized electrons are ne- 
glected, (ii) Because the Zn-3d electrons are localized, 
they are not well described within the local spin den- 
sity approximation (LSDA) to density functional theory 
(DFT). To treat these electronic states adequately, one 
has to go beyond the LSDA, for example by applying 
the self-interaction correction (SIC) to the LSDA. (iii) 
Previous studies often focused on the WZ and RS struc- 
tures of the ordered alloys (i. e. a; = or 1). A detailed 
investigation of the complete transition path (with con- 
tinuous variation of x) is still missing. For the latter it is 
inevitable to treat disordered alloys, for example within 
the coherent potential approximation. 

In the present paper, we report on a systematic first- 
principles study of structural and electronic proper- 
ties of ordered and disordered Mga;Zni_^0 alloys us- 
ing an all-electron full-charge density Korringa-Kohn- 
Rostoker (KKR) method4^ii^ Alloying of MgO and ZnO 
is described within CPA^ as formulated in multiple- 
scattering theory (KKR-CPA)4^ One objective of this 
investigation is to describe accurately the structural 
phase transition from the WZ to the RS structure in 



Mg3;Zni_^0, by continuously increasing the Mg concen- 
tration. The second objective addresses the formation of 
the fundamental band gap in Mga;Zni_a;0. It is further 
shown how the hybridization between Zn-Sd and 0-2p 
states affects the width of the band gap, by comparing 
results obtained within the LSDA with those obtained by 
applying the self-interaction correction4^i^ The validity 
of the present approach, especially the use of the CPA, is 
discussed by comparing our results with those of previ- 
ous studies i^i^i^ In summary, our study addresses im- 
portant issues in Mg^^Zni-^jO alloys which, with respect 
to magnetic tunnel junctions, might also be relevant for 
magnetoelectronics. 

The paper is organized as follows. The model of 
the transition path from the WZ to the RS structure 
is sketched in Section [III Details of the computational 
approach are presented in Section IIIIl Results are dis- 
cussed in Section IIVI By comparing our results for the 
ordered alloys with those of other theoretical work and 
with experiment, the validity of our approach has been 
established. Our main results for the disordered alloys 
are discussed in Section |Vl In Section IVll the formation 
and evolution of the fundamental band gap as a func- 
tion of Mg and Zn contents in the system is analyzed. 
Concluding remarks close the paper in Section IVlII 



II. MODELING THE STRUCTURAL WZ ^ RS 
PHASE TRANSITION 

To investigate structural phase transitions in 
Mga;Zni_i^O, a transition path suggested for the 
continuous structural WZ <-> RS transformation in GaN 
was adoptedi^ This path has been already successfully 
applied in studies of structural deformations in MgO^ 
and in ZnO 4^ The structural transition is described 
as a homogenous strain deformation from the WZ to 
the RS phase by passing an intermediate hexagonal 
structure (HX, refered to as h-MgO in Ref. [HI). This 
hexagonal structure can as well occur in epitaxial 
systems due to the reduction of the interlayer distance, 
which was recently observed for thin ZnO filmsj ^^i^° 
A similar scheme is the Bain's path for the transition 
from the face-centered-cubic to the body-centcred-cubic 
structurci^ 

In the first step on the transition path (WZ ^ HX), 
the internal parameter u is linearly increased while simul- 
taneously decreasing the '^/a ratio (Fig.[T]). At u = 1/2 the 
space group changes from PG^mc to PG^/mmc. In the 
second step (HX ^ RS), the lattice is compressed uni- 
axially along the [1010] direction and simultaneously '^/a 
is decreased further. The WZ phase is characterized by 
a = b, c/a = -y/^, V = 1/3 and u = s/s, whereas the HX 
phase has a ~ b, '^■/a — 1.2, v = 1/3, and u = 1/2. For the 
RS phase, a = 6 = c and w = ?i = 1/2. For details, see 
Ref. [H. 

To perform total-energy minimizations (see Sec- 
tion [lTl|, the appropriate unit cell of an orthorhombic 
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FIG. 1: (Color) Structural phases in Mg^Zni_j;0. The 
wurtzite structure (a: WZ, left column), the intermediate 
hexagonal structure (b: HX, central column), and the rock- 
salt structure (c: RS; right column) are shown in top view 
(top row) and in side view (central row). The hexahedral 
chemical units are depicted in the bottom row. The lattice 
parameters a, b, and c for the common space group Cmc2i axe 
indicated. The ly-plane of unit cells for PGamc, PQz / mmc, 
and FmSm space groups for WZ, HX, and RS respectively are 
displayed in the top row in magenta, u describes the internal 
displacement of cations sites (red, dark spheres) with respect 
to the central planes of the polyhedra. For u = 1/2, metal 
sites lie within the central planes which are spanned by the 
oxygen sites (yellow, bright spheres). 



lattice with space group Cmc2i was used (Fig. [T|). The 
latter is a common subgroup of all three phasesi^ In our 
calculations, the parameter v, specifying the relative po- 
sitions of the subsublattices^^ and the ratio b/a (defined 
in Fig.[T]) are fixed by the geometry of a particular phase. 
The parameter u, which is the internal vertical displace- 
ment of the cation atoms in the oxygen plane, and the 
c/a ratio were obtained by the total energy minimization 
for the pure compounds (i.e., WZ-ZnO and RS-MgO) 
within LSDA. They are kept fixed at these values for 
Mg3;Zni_2:0 alloys for all concentrations x. 



III. COMPUTATIONAL DETAILS 

The electronic and geometric structures of 
Mga;Zni_2^0 alloys arc obtained within density- 
functional theory. The local density approximation 
and its self-interaction correction (see below) are 
implemented in a multiple-scattering approach (the 
Korringa-Kohn-Rostoker method, KKR)iS Disordered 
alloys are described within the coherent potential 
approximation (CPA)i^ 

For closed-packed systems, for example transition met- 
als, the crystal potential is commonly approximated as a 



sum over the so-called 'mufSn-tin potentials— (that is, 
spherically symmetric potentials centered at each lattice 
site, while in the interstitial region the potential is con- 
stant). For open systems, this approximation results in 
a relatively poor description of the electronic structure. 
It has turned out that the usual trick of inserting the so- 
called 'empty spheres' into the interstitial region is not 
sufficient for MgO and ZnO, in particular in view of the 
high degree of accuracy needed in the evaluation of total 
energies. Obviously, the latter is inevitable for a reli- 
able description of structural phase transitions. Conse- 
quently, a full-charge density approximation was applied. 
Here, the total energy is estimated from the non-spherical 
charge density and the full potential. All radial integrals 
are calculated using the unit cell geometry; the single-site 
problem is still solved with a spherical potential averaged 
within a Voronoi cell. The accuracy of this approach is 
as good as that of a full-potential method but is not as 
time-consuming. The validity of the full charge-density 
approximation was checked by comparison with results 
of the corresponding KKR full-potential calculations for 
the ordered WZ-ZnO and RS-MgO compounds. 

The disordered Mg3;Zni_a;0 can be viewed as a substi- 
tutional binary alloy in which the metal sublattice (Zn, 
Mg) is subject to chemical disorder. The oxygen sublat- 
tice remains unaffected. Hence, the coherent potential 
approximation (CPA) is an obvious choice for describing 
M gxZni_2,0 at an arbitrary concentration x. Within the 
CPA, Mg and Zn impurities are embedded into an effec- 
tive coherent potential medium which is determined self- 
consistentlyi^i^ Since the KKR-CPA is a single-site ap- 
proximation, it does not allow to investigate the influence 
of chemical short-range order on the electronic properties 
(SRO can be taken into account by the cluster CPA;^ 
the locally self-consistent Green function method^ or the 
non-local CPA approach^). However, for many alloys 
the single-site CPA provides a reasonable description of 
the electronic structure.— For Mga;Zni_2,0, the valid- 
ity of the KKR-CPA was established through a direct 
comparison with KKR supercell calculations (which in- 
clude SRO effects) for x = 0.25, 0.50, and 0.75. The 
major structural and electronic properties, such as the 
equilibrium lattice constants, bulk moduli, equilibrium 
pressures, and fundamental band gaps, are reproducible 
within both the KKR-CPA and the supercell approach. 
Hence we conclude that SRO effects cannot be ruled out 
but they are of minor importance for the issues addressed 
in this work. A detailed investigation of SRO effects is 
beyond the scope of the present work. 

Although the Mg-2p electron states lie comparably 
deeply in energy (about 2 Rydberg below the valence 
bands), they have been treated as valence states. It is 
found that the hybridization of the associated electronic 
states with the valence states is important for an accu- 
rate determination of total energies, especially for the 
evaluation of both the °/a ratio and the parameter u in 
the WZ structure. The same procedure was applied in an 
earlier all-electron study of MgO, using a full-potential 
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linearized mufRn-tin orbital (FP-LMTO) methodi^ 

Since the Zn-Sd states are energetically close to the 
0-2p states, they have to be considered as valence elec- 
trons as well. It is well-known that the hybridization 
of the localized Zn-Sd electrons with the 0-2p electrons 
is crucial for an accurate investigation of the band-gap 
formation in ZnOi^ Being strongly localized, the Zn-Sd 
states are not adequately described within the local den- 
sity approximation. The LSDA contains the (unphysical) 
self- interaction of an electron with itself (see e. g. Refs.l47l 
and IgU ) , an effect whose importance increases with the 
degree of electron localization. As a consequence of the 
self-interaction, the energy levels of the localized elec- 
trons obtained within the LSDA lie too high and their 
hybridization with the 0-2p states is too strong. The 
calculated fundamental band gap in ZnO is therefore con- 
siderably too small as compared to experiment. 

An improved description of localized electrons is 
achieved by the self-interaction correctional^ in which 
the unphysical self-interaction is removed from the LSDA 
exchange-correlation functional. In this approach, SIC 
is applied to various configurations of localized electron 
states, and the configuration with the lowest total energy 
defines the ground state energy and configuration. The 
SIC-LSD approach, being based on a variational princi- 
ple, is parameter free and treats on equal footing both 
itinerant and localized electrons. When no localized elec- 
trons are present in the system, then the SIC-LSDA to- 
tal energy is equivalent to the LSDA total energy. Thus 
the LSDA energy functional is a local minimum of the 
SIC-LSDA functional. For ZnO in both the RS and the 
WZ phase it is found that all 3d-electrons have to be 
Sl-corrected. Since the Zn-3d states have semi-core char- 
acter, application of the SIC leads to a uniform increase 
of the binding energies of these states, to a decrease of 
the total energy, and to an increased fundamental band 
gap in ZnO. Because other electronic states were treated 
within the LSDA, we label this calculations SIC-LSDA 
for short. 

To check whether SIC affects only the fundamental 
band gap or also structural properties, we have compared 
results of a SIC-LSDA calculation for ZnO with those of 
an LSDA calculation, based on the same computer code. 
It has been shown that structural properties calculated 
within both LSDA and SIC-LSDA are very close, whereas 
the band gap obtained within SIC-LSDA has been sub- 
stantially larger than in the LSDA case. This finding has 
been further verified for disordered Mga;Zni_j;0 alloys at 
selected concentrations. As a result, in order to save com- 
putational costs, only LSDA has been used for the studies 
of the structural properties. In contrast, here, the funda- 
mental band gap has been investigated at the respective 
equilibrium lattice constants within SIC-LSDA. For both 
the LSDA and SIC-LSDA calculations the Perdew-Wang 
exchange-correlation functional has been appliedi^ 

The equilibrium volumes, bulk moduli, pressures, and 
enthalpies have been calculated for zero temperature 
from the total energy fitted to the Murnaghan equa- 



tion of stateiSs Lattice vibrations, finite-temperature ef- 
fects, and relativistic corrections have not been consid- 
ered. The angular momentum cut-offs ^max have been 
chosen as 3 for the Green function expansion and 6 for 
both the charge-density and the potential representation. 
The convergence with Zmax and with normalization of 
the Green function was significantly improved by use of 
Lloyd's formula in this work i^^i^'^ The concentration, x, 
of Mg impurities in ZnO has been varied in steps of 5 %. 



IV. STRUCTURAL PROPERTIES OF 
ORDERED ZnO AND MgO IN RS AND WZ 
PHASES 

Due to the large number of structural parameters, a 
complete structure optimization of Mga;Zni_2;0 alloys is 
an involved task. Apart from these parameters, the alloy 
composition x is an additional degree of freedom which 
complicates the problem further. Therefore, the number 
of parameters to be optimized has been reduced by con- 
centrating on volume changes in the WZ, HX, and RS 
structures upon variation of x (see Section HI)) . 

To our knowledge, the KKR method was not used be- 
fore for optimizing WZ and HX structures. Hence ex- 
tensive calculations for pure ZnO and MgO have been 
required to determine both the optimum u and ^/a. In 
accomplishing this, the procedure suggested in Ref. [s^ 
has been followed. Since the calculation of lattice relax- 
ations from forces is rather complicated and not suffi- 
ciently accurate within the KKR method, we have cal- 
culated the total energies consecutively varying three pa- 
rameters: the lattice constants a, the '^/a and the internal 
parameter u. The total energy has been calculated at a 
given u for a sufficiently dense mesh of lattice constants 
and then fitted to the Murnaghan equation of state. The 
initial value of the parameter u has been taken from the 
experiment for a particular structure and then succes- 
sively varied as long as the absolute total energy min- 
imum has not been reached. To establish the validity 
of the present KKR approach for structure optimization, 
our results have been compared with those of other first- 
principles calculations. We note that our calculations 
for both ZnO and MgO in various structures have been 
carried out for the same unit cell on the same level of 
approximation. 

The adequacy of the present approach is evidenced by 
the agreement with experimental data and the results of 
other theoretical approaches (see Tables U and [n| . It is 
found that our results agree with those obtained by the 
other all-electron methods, as opposed to the pseudo- 
potential methods. In particular, the parameters ob- 
tained in this work compare well with those reported by 
Limpijumnong and coworkers i^^ii^ This agreement might 
be attributed to the fact that both studies follow the 
same optimization scheme (as suggested in Refs. [s^ and 
lish and that the Zn-3d as well as the Mg-2p electrons 
have been treated as valence electrons. 
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TABLE I: Properties of ZnO in the wurtzite (WZ), tire hexagonal (HX), and the rock-salt (RS) phase. Equilibrium volumes 
Vo, '^ja ratios, internal parameters ti, and bulk moduli _Bo are compiled for theoretical (present work displayed in bold) and for 
experimental works. PPW, LCAO, and FP-LMTO are short for pseudo-potential plane waves, linear combination of atomic 
orbitals, and full-potential linearized mufBn-tin orbital, respectively. 



Phase 




Theory 


Experiment 


WZ 


V, (A^) 


22.80", 22.83, 22.87', 22.91", 22. 93% 23.4"-', 23.62% 23.78" 


23.80 "' 






1.590% 1.602, 1.605", 1.607", 1.608^, 1.610", 1.614% 1.617^* 


1.602'*' 




U 


0.379% 0.380"s% 0.381" 


0.382''' 




Bo(GPa) 


154, 154", 155", 157-', 160% 162"'"' 


143'', 183' 


HX 


Vo (A^) 


22.12 








1.200 






Bo(GPa) 


165 




RS 




18.70", 18.76", 18.87% 18.88, 18.98% 19.08% 19.45" 


19.48', 19.60" 




Bo(GPa) 


200", 201, 203", 206% 210", 211% 219» 


202', 228'' 



"Ref. 


4g 


PPW. 


''Ref. 


3C 


LCAO. 


'^Ref. 


36 


LCAO. 


"^Ref. 


41 


PPW. 


"^Ref. 


4C 


PPW. 


/Ref. 


66 


FP-LMTO. 


fRef. 


67 


FP-LMTO. 


"Ref. 


27 




'Ref. 


H 





TABLE II: As Table H] but for MgO. Missing references are given in Table H] 



Phase 




Theory 


Experiment 


WZ 


w 

Bo(GPa) 


22. 50^,, 22.53", 23.15% 23.20 
1.550", 1.610, 1.620% 1.633-' 
0.380, 0.380 ^ 
121, 131", 137' 




HX 


Vo (A^) 
Bo(GPa) 


20. 90^, 21.71 
1.200, 1.200' 
135, 148' 




RS 


V, (A^) 


17.54", 17.80% 17.96% 18.03% 18.19, 18.65"^ 


18.67', 18.75"' 




Bo(GPa) 


167, 170", 172% 174^*, 178% 186'' 


169", 172° 



mef. 
*^Ref. 
'Ref. 
™Ref. 
"Ref. 
°Ref. 



FP-LMTO. 
FP-KKR. 



As for the structural properties, we find agreement con- 
cerning the volume ratios of the different phases and the 
equilibrium pressures at the phase transitions (see Ta- 
bles IIIl and |TV| . For ZnO there is only the WZ ^ RS 
phase transition while for MgO there are two transitions, 
namely RS <-> HX and RS <-> WZ. In particular for ZnO, 
our theoretical results compare well with experimental 
ones (see Table Hill . For MgO we find agreement with 
other theoretical estimations too (Table [IV|. However, 
the equilibrium pressure for the RS-HX phase transition 
might be an exception. The apparent disagreement with 
the result reported in Ref. [s^ may be related to the insta- 
bilities of the WZ and the HX structures (Figs.[2]and[3^). 
Moreover, due to the negative pressure these phases can- 
not be realized experimentally, hence ruling out a clar- 



TABLE TIL Structural phase transition in ZnO. Given are the 
volume ratio and the equilibrium pressure at the WZ ^ RS 
phase transition. References as in Table |T] Present work 
displayed in bold. 





Theory 


Experiment 


Vwz/Vrs 


1.20% 1.21, 1.22""", 1.249 


1.21', 1.22" 


Pwz^RS (GPa) 


3.9", 6.6", 8.0% 8.2", 8.6, i. 


.7" 8.7", 9.1' 



ifying comparison of theory with experiment. In favor 
of our work, we would like to mention that the present 
calculations reproduce well '^ja and u. One might spec- 
ulate that earlier implementations of multiple-scattering 
theory failed to optimize these structural parameters in 
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TABLE IV: As Table HITl but for MgO. For MgO there are 
two transitions, RS <^ HX and RS ^ WZ. Experimental data 
are not available. References as in Table |ll] Present work 
displayed in bold. 



Theory 


Hhx/Vrs 


1.17^ 1.19 


Prs-hx (GPa) 


-16.2^ -8.5 


Vwz/Vrs 


1.24, 1.26^ 1.28" 


Prs^wz (GPa) 


-11.1, -8.4^ 



any open structure due to the slow angular momentum 
convergence of the Green function (see Section IIIIl the 
convergence is significantly improved by using Lloyd's 
formula) 

In summary, the agreement of our results for the or- 
dered alloys with those of other theoretical investigations 
and experiments indicates that the present approach is as 
reliable and accurate as any other state-of-the-art meth- 
ods. 



V. STRUCTURAL PROPERTIES OF 
DISORDERED Mg^Zm^O ALLOYS 

Having established that our results for ordered ZnO 
and MgO are consistent with the structural properties 
obtained by other theoretical methods and by experi- 
ments, we now turn to the discussion of results for dis- 
ordered Mga,Zni„j.O alloys. To reduce the number of 
parameters to be optimized, both ^/a and u were fixed 
respectively to 1.6 and 0.38 in the WZ structure and to 
1.2 and 0.5 in the HX structure. This approximation is 
justified by the weak dependence of both '^/a and u on the 
atomic species (Mg or Zn; cf. Tables U and |lT| . Conse- 
quently, the structure (the actual phase), Mg concentra- 
tion X, and unit cell volume remain to be varied. The 
results of the total-energy calculations involving these 
three variables are discussed below. The results of the 
total energy calculations are comprised in the formation 
enthalpy (at T=0) 

Aiy„(Mg,Zni_,0) = i;„(Mg,Zni_,0) 

-x£;Rs(MgO) - (1 - x)£;wz(ZnO) ^ ' 

of a structure a of Mgj;Zni_a;0, relative to the most sta- 
ble forms of ZnO (WZ) and MgO (RS) compounds,^ 

The formation enthalpy is positive throughout, as 
seen in Fig. [21 in agreement with previous theoretical 
studiesj^i^ There are two global minima (i. c. the most 
stable phases; black, dark regions), for pure ZnO in WZ 
{x = 0) and pure MgO in RS phase {x = 1). This finding 
implies a tendency towards phase separation if the inte- 
gration of different constituents into the medium cannot 
be maintained)^ 

In the following, five special cases are discussed in more 
detail, namely x = 0.0, 0.33, 0.52, 0.71, and 1.0. 



Enthalpy per formula unit (eV) 




16 18 20 22 24 26 28 

Volume per formula unit (A^) 



FIG. 2: (Color) Formation enthalpy of Mgj;Zni_j;0 alloys 
versus equilibrium volume and Zn concentration, depicted as 
color scale [cf. eq. ([T|]. Straight bold lines mark separations 
between the phases (WZ, HX, and RS). 



With low Mg concentration, the WZ phase is favorable 
in any case (Fig. [3^). The pressure needed for a WZ <-> 
RS transition is positive and increases monotonously 
(Fig.©. 

For X < 0.33 the total energy of the HX structure is 
larger than for the RS and the WZ phases, implying that 
an intermediate HX phase cannot be established (Fig. [2]). 
Consequently, a direct WZ ^ RS transition is possible at 
positive pressure. At a; = 0.33 (Figs. [HJa and \^ this 
phase transition can take place at zero pressure. This 
finding is consistent with the theoretical work of Sanati 
and coworkers^l and is also observed experimentally in 
ZnO-MgO heterostructuresji 

At X = 0.71 the WZ <-> HX and the HX ^ RS transi- 
tions occur at the same pressure (Fig.|3]). This is evident 
from Fig. [5}; because all total energy curves can be con- 
nected by a single tangent. We note that for x = 0.52 
the WZ <^ HX transition is possible at zero pressure, and 
the corresponding total energies relative to RS phase are 
identical (see Fig. [5]). 

For pure MgO {x = 1.0; Fig. [3ji), the RS structure 
exhibits the global minimum at a volume of 18.19 
which is about 2 % less than the experimental value (see 
Table [lT|. By applying a negative pressure (i.e. by in- 
creasing the volume along the tangent) the RS phase is 
transformed into the HX structure, in agreement with the 
work of Limpijumnong and coworkers 

According to the phase diagram (Fig. [2]) , the HX phase 
is an intermediate phase between the WZ and RS struc- 
ture for 0.33 < a; < 1.00. Although the HX ^ RS transi- 
tion pressure increases with Mg content, it remains neg- 
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Volume per formula unit V (A^) 



FIG. 3: (Color online) Relative total energies of Mg^Zni_^0 
alloys in RS (black solid line), HX (green dashed line), and WZ 
(blue dotted line) structures at four different concentrations: 
(a) X = 0.0 (pure ZnO), (b) x = 0.33, (c) x = 0.71, and (d) 
X = 1.0 (pure MgO). Total energies in cases (a) are taken 
relative to the WZ phase at a; = 0.0 and, in (b), (c), and 
(d), relative to the Mg phase at a; = 1.0. Red arrows show 
experimental values of volume (see Tables HI and 



FIG. 5: (Color online) Total energies of Mga;Zni_a;0 alloys 
in the HX (blue dashed line) and the WZ (green dotted line) 
structures, taken relative to the total energy of the RS struc- 
ture at equilibrium volumes. 



ative in the whole range of concentrations (Fig. 2]) . 



VI. BAND GAP IN Mg^Zni_^0 ALLOYS 
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FIG. 4: (Color online) Equilibrium pressure in Mg^Zni-^O 
alloys: RS <-» WZ (green dashed line) and RS <-> HX (blue 
dotted line). 



Since Mga;Zni_a;0 alloys arouse great interest as band 
gap-engineering materials, in this section we study the 
evolution of their fundamental band gap as a function of 
the Zn and Mg contents. We restrict our consideration to 
the RS and WZ phases, as only these appear interesting 
from the experimental point of view. 

Band gap formation can be correctly described only 
within a many-body theory which takes properly into 
account the electron-hole excitations. However, the 
GW approximation^, probably the most popular first- 
principles many-body approach for calculating excitation 
energies, is too expensive for studying alloys with arbi- 
trary concentrations and also having constituents with lo- 
calized electrons. Therefore, for this study we have used 
only the LSDA and SIC-LSDA approaches. The calcu- 
lated band gaps have been estimated from the density of 
states (DOS), band structure and Bloch spectral func- 
tion, corresponding to the theoretical equilibrium lattice 
structure of a given chemical composition. 

Using LSDA for RS-MgO we obtain a band gap of 
5.15 eV, which is 67 % of the experimental value (7.7 eV). 
For the pure WZ-ZnO, LSDA gives the band gap of 
0.8 eV, which is only 23% of the experimental value 
of 3.35 cV. In MgO, the difference between the calcu- 
lated and experimental values originates mostly from the 
fact that the electron-hole excitations arc neglected in 
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LSDA, which can be corrected effectively with the GW 
approximation^^ For ZnO the problem is more com- 
plex, as in addition to the neglect of the electron-hole 
excitations, the localized nature of Zn 3d electrons is 
not adequately represented within LSDA. This failure of 
LSDA is largerly related to the inherent unphysical self- 
interaction. The latter affects the band gap by placing 
the localized 3d electrons of Zn at too low binding ener- 
gies, thus leading to their strong hybridization with the 
0-2p states. In LSDA the Zn-3d states are about 3eV 
too high with respect to the experimental value of about 
7.5 eV-8.5 eV^^ (see the DOS of ZnO in Fig.©. In the 
RS and WZ phases, the 3d states are placed at slightly 
different energies, and in the RS the corresponding band 
width is larger due to the more closed-packed crystal en- 
vironment. The strong pd- hybridization dramatically re- 
duces the band gap, which can be improved only slightly 
within the GW approximation, if based on the LSDA 
Green's functioni^ It is the aim of our future studies to 
use the SIC-LSDA band structure of ZnO for the sub- 
sequent GW calculation. The SIC-LSDA band gap we 
have calculated is 2.42 cV for WZ-ZnO and 2.96 eV for 
RS-ZnO. Thus, the SIC-LSDA band gap for the WZ-ZnO 
constitutes 69% of the experimental value of 3.35 eV. 
This tells us that SIC-LSDA is at least as good for the 
pure ZnO, as LSDA is for the pure MgO. The reason be- 
ing, that SIC-LSDA describes both itinerant and local- 
ized electrons on equal footing. Of course, in MgO the 
band gap is constituted by the sp electrons which are itin- 
erant and thus unaffected by the spurious self-interaction, 
while in ZnO the Sd electrons of Zn play a defining role 
in establishing the band gap. The SIC-LSDA approach, 
by removing the unphysical self-interaction of all 10 Zn 
d electrons, describes them much more adequately than 
LSDA. What happens as a result of SIC is that all the 
d-states of Zn move to higher binding energies and the 
Zn-derived bands become narrower thereby reducing the 
pd-hybridization and increasing the band gap, as seen 
in Fig. [51 The structural properties of Mga;Zni_j;0 are 
believed to be little affected by this uniform shift down- 
wards in energy of the d-states, leading mostly to a uni- 
form lowering of the total energy. 

Being an effective one-electron ground state theory, 
SIC-LSDA does not provide a quasi-particle spectrum to 
compare with spectroscopies. Missing the crucial screen- 
ing/relaxation effects (self-energy), it predicts the 3c? Zn 
electrons at too high binding energies, as opposed to 
LSDA where they come out too low (Fig. [6]). One can 
implement a simple fix to correct for the screening ef- 
fects in SIC-LSDA, based on the Slater's transition state 
theoryi^SiiS Following the concept of the latter, we cal- 
culate the SIC-LSDA-based removal energies of localized 
electrons as an average of the calculated SIC-LSDA and 
LSDA d-state expectation values 

stsa = AHlsda + Vsic\d > + < d\HLSDA\d >) . 

(2) 



1 ' I ' I ' I ' I ' r 




Energy (eV) 

FIG. 6: (Color) Density of states (DOS) of ZnO in the RS (a) 
and the WZ (b) phase obtained with the LSDA (black solid 
lines), the SIC-LSDA (red dashed lines) and SIC-TSA (blue 
dash-dotted lines) aooroaches. Insets show the behaviour of 
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0.0 0.2 0.4 0.6 0.8 1.0 
Concentration of Mg 

FIG. 7: (Color online) Band gaps of MgjjZni-^O alloys in 
the RS and WZ phases versus Mg content, calculated within 
LSDA and SIC-LSDA (line styles as indicated). 



Effectively, the above equation states that only half of 
the SIC potential should be applied at the stage of cal- 
culating the density of states, after the self-consistency 
has been achieved. We refer to Eq. [2] as the transition 
state approximation (TSA) and show the resulting DOS 
of ZnO in Fig. [HI We can see that Zn d-states, calculated 
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TABLE V: The band gap Eg for MgO (RS), ZnO(WZ) and 
the bowing parameter for Mga;Zni_a;0 alloys. 





LSDA 


TSA SIC-LSDA 


Experiment 


Band gap [eV] 


RS-MgO 


5.15 




7.70 


WZ-ZnO 


0.80 


1.74 2.42 


3.35 


Bowing parameter [eV] 


RS 


1.90 


0.91 0.59 


0.70± o.ais 


WZ 


4.88 


4.05 3.65 


3.60± 0.6^ 



using TSA, appear at lower binding energies, as com- 
pared to the strict SIC-LSDA result. Consequently, the 
TSA-binding energies of Zn d states are in better agree- 
ment with experimental values of 7.5 eV-8.5 eV— i^. The 
effective hybridization of Zn d-states with the oxygen p- 
states is stronger as in the SIC-LSDA case, which leads to 
a reduction of the band gap to 2.32 eV and 1.74 eV (52% 
of the experimental value) in RS and WZ phases, respec- 
tively. 

In Mga;Zni_j;0 alloys the size of the band gap depends 
on the concentration of Mg impurities, as seen in Fig. [7] 
and Table El While at the Mg rich end the LSDA band 
gaps are closer to experiment, at the other end the SIC- 
LSDA band gaps are more adequate. As for the TSA 
results, they fall mostly in between the LSDA and SIC- 
LSDA band gaps, especially for small Mg concentrations, 
where the hybridization of the O 2p bands with the Zn 
3d states is of great significance (see Fig. [T]). The be- 
haviour of the band gap as a function of Mg concentration 
is also different between the various approaches. While 
the LSDA curves are rather parabolic, the SIC-LSDA 
and TSA band gaps seem to change almost linearly with 
concentration, thus following more closely experimental 
resultsiSiii^. This is mostly due to the fact that the actual 
magnitudes of the band gaps change slower with concen- 
tration in SIC-LSDA and TSA than in LSDA. The reason 
being that the larger Zn content, the more inadequate 
LSDA is and the smaller the resulting band gaps. 

Although, both LSDA and SIC-LSDA underestimate 
the size of the band gap. Eg, its dependence on the Mg 
concentration, x, can be compared with experimental re- 
sults. One way of doing it is to estimate the so-called 
bowing parameter, b, appearing in the commonly used 
definition of the fundamental band gap dependence on 
the composition x, namely 

Eg{x) = a;£:g(MgO) + (1 - x)Eg{ZnO) - bx{l - x), 

where the bowing parameter is given by 

b = 2i;g(MgO) + 2Eg{ZnO) - 4Eg{Zno.5oMg^,,^0). (3) 

Thus to evaluate the bowing parameter one needs to 
know the band gaps of the pure MgO(RS) and ZnO(WZ), 
as well as their alloy with the concentration x — 0.5. 

According to our total energy calculations the 
Mga;Zni_2:0 at X = 0.5 occurs in the RS phase which 



is in agreement with the experiment of Chen et al.^ 
However, the crystal structure of Mga;Zni_a;0 thin film 
alloys depends strongly on the growth method. In vari- 
ance to Ref. [3 the Mg^Zni_a;0 films with x < 0.53, 
prepared with the PLD procedure, were found to have 
the WZ structure. Therefore, for direct comparison with 
the experiments we have used Mgo.50Zno.50O in both RS 
and WZ phases. The resulting calculated bowing param- 
eters for Mg^Zni_2:0 alloys are presented in Table IVl We 
can see that the LSDA systematically overestimates the 
bowing parameter for both structures, while the straight 
SIC-LSDA agrees very well with experiments. The TSA 
results fall in between the LSDA and SIC-LSDA values. 

Summarizing this section, we have to say that despite 
rather good agreement of our SIC-LSD bowing parameter 
with experiments, to be truly predictive, one would need 
a more robust method like a combination of SIC and GW. 
This way we could also predict the correct magnitudes of 
the band gaps and enter the serious business of band gap 
engineering. 

VII. CONCLUSIONS AND OUTLOOK 

Structural phase transitions and the fundamental 
band gaps of Mga;Zni_a;0 alloys have been investigated 
by detailed first-principles calculations. The multiple- 
scattering theoretical approach used here (Korringa- 
Kohn-Rostoker method) allows to treat disordered al- 
loys within the coherent potential approximation, that 
is Mga;Zni_a;0 alloys with arbitrary Zn concentration 
X, with an accuracy as good as in other first-principles 
methods. The importance of treating localized states of 
ZnO appropriately within the framework of the local- 
density approximation to density-functional theory is es- 
tablished, thereby confirming the usefulness of the self- 
interaction correction. 

The delicate interplay of geometry and electronic 
structure is not only of importance for bulk systems, as 
shown in this work. In nanotechnology, interfaces and 
surfaces play an essential role. Therefore, a correct de- 
scription of disordered alloys and their electronic struc- 
ture, in particular at the Fermi energy, is necessary for 
predicting material properties, besides confirming and 
explaining experimental results. However, to be fully 
predictive regarding the properties of such systems as 
M ga;Zni_3:0 alloys, of importance for device applications, 
one needs a first-principles approach like a combination 
of GW with SIC-LSDA. 
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